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^D . A detailed calculation of the finite temperature dynamic structure function 
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of the free Bose gas is presented and discussed. After a short derivation of the 



B 

i-rt ■ expressions describing the exact response S{q,uj;T) and the Impulse Approx- 

O ■ imation (lA) SfAiQy '-^j t)i their main properties and their evolution with q and 



T are analyzed. The lowest order energy weighted sum rules of both S{q, uj; T) 
and Sia{(1t^'iT) are also derived and commented. Finally, the g — > oo asymp- 
totic behavior is analyzed and discussed in terms of scaling laws. 
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Neutron scattering on systems like pure "^He has proved to be among the most fruit- 
ful methods to gather relevant information about the fundamental properties of Bose con- 
densed systemslil. Although high momentum transfer neutron scattering is known to provide 
much information about the structure of the Bose condensate, no particular formalism has 
been able to accurately describe how temperature influences the dynamic structure func- 
tion S{q, uj) and how the temperature dependence of the condensation fraction value affects 
S{q, u) at high and low momentum transfers. More recently, the experimental discovery of 
Bose condensation in clouds of alkali atoms has revived the interest in the analysis of free 
Bose systems and the way in which Bose-Einstein condensation sets inBu. In previous works 
we presented detailed calculations of the T = and T > dynamic structure function of the 
free Fermi gas, their coherent and incoherent parts and the way in which they evolve with 
momentum transferQ'S. In this work we extend this analysis to free Bose systems, focussing 
on the evolution of S{q, u) and its coherent and incoherent parts with q and T, the role of 
the Bose-Einstein condensate in the response and finally on scaling laws. 

The dynamic structure function S{q, u; T) describes the way in which a system of A^ 
particles responds to a probe which produces a density fluctuation modelized through the 
operator pq = I]jLi e*'''"^ . at finite temperature this is given by&El 

S{q,uj-T)= Y: ^e-^(^"-^^)l|(m|p, |n)r<5(E^-E.-u;) , (1) 

{n,m} 

where q and u are the momentum and energy transferred by the probe to the system. 
Z = J2{n} e~^^^"~^'^^ stands for the Grand Canonical partition function, while /5 = l/fc^T 
is the inverse of the temperature and fi is the chemical potential. 

The dynamic structure function is related by Fourier transformation to the density- 
density response function S{q, t) which may in turn be separated in its coherent and inco- 
herent parts 
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H and r^ being the Hamiltonian of the system and the position operator of particle j, respec- 
tively. With these definitions the incoherent response takes into account only the diagonal 
terms of the density-density response function, while the coherent response incorporates the 
contribution coming from all the other terms in S{q,t). Notice that these definitions are 
not unique as in fact frequently Scoh{(l,t) is identified with S{q,t) while the last function in 
Eq. (^ becomes simply the difference between the new Scohi^yt) and Sinc{q,t). 

The dynamic structure function is also related to the imaginary part of the dynamic 
susceptibility x{Qi^'-,T). Taking into account that x{q-,^\T) = — 7r(l — e^^'^)S'(g, cu; T)ffl, 
one finds 
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that is written in terms of the dimensionless variables q = q/ko, u = w/eo and T = T/eo, 
where ko = p^/^ and cq = /i;g/2m define the momentum and energy scales. Finally, in Eq. (^ 
z stands for the fugacity that is related to the dimensionless chemical potential through 
z = e^/^, where jl = /i/eo. This last quantity can be derived recalling that at equilibrium 
particle number conservation imposes 



z~^e^ — 1 
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where no(T) is the condensate fraction value at temperature Tb. For the free Bose gas, at 



T <Tr 



no(T) = 1 






3/2 



(6) 



where Tc ~ 6.63 is the value of the dimensionless Bose-Einstein transition temperature 
which in the current dimensionless scheme becomes density and mass independent. 



Notice that the first term in Eq. (^ is proportional to the condensate fraction value no{T) 
and thus characteristic of Bose liquids: when the system is given a net momentum transfer q, 
all particles in the condensate respond equally to the perturbation and the response peaks at 
a single frequency u = (f. On the other hand, the second term is the contribution of all the 
other particles not laying in the condensate and is similar to the partticle-hole contribution 
to the total response of Fermi systems. 

In the free case and due to the absence of an interatomic potential, the incoherent 
response coincides with the Impulse Approximation, 



Sinc{q,uj]T) = SiA{q,uj]T) = j^^ J dkn{k) 6(^(k + qf -V^ -iS). 

Using the momentum distribution of the free Bose gas one gets the following result 
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Sinc{q,Co;T) = no{T)6(uj-q^^ 
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while Scoh{Q,^]T) is the difference between (^) and (||). 

Below the Bose-Einstein transition temperature, where z = 1, one can introduce a 
new set of variables Q = q/T^/'^ and u = u/T, such that the ratio of the non-condensate 
contributions and T^/^ becomes temperature independent for all three responses when Q 
and u are taken as new independent variables. Therefore, it is useful to introduce three new 
structure functions that take this fact into account, 
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ScohiQ, y) = S{Q, U) - Sinc{Q, V 
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where the subscript nc indicates that only the non-condensate contribution to the different 
responses has been taken into account. 

The total, coherent and incoherent non-condensate responses at T < T^ are plotted in 
Fig. (j^) for Q = 0.7, Q = 1.5 and Q = 2. Both the q dependence and the T dependence of the 



responses are represented in there, because Q can be increased either by increasing q at fixed 
T or by lowering T at fixed q. In the first case, the responses at different Q can be directly 
compared, while in the latter case one should bear in mind that Sndq, ^] T) = T^/'^S{Q, v). 

As it is readily seen from Eq. (^, both the total and the coherent responses of the free 
Bose gas presents logarithmic singularities aX u = if and Cj = —f, (that is, sX u = iQ^) 
when T is lower than the Bose-Einstein transition temperature, as in this case z = 1. On 
the other hand, at these temperatures SmdQ, i^) presents only one singularity located at 
V = Q^. In this way, the pole aX v = Q^ has contributions coming from both the coherent 
and the incoherent responses, while the pole at v = —Q^ is of a completely coherent nature. 
As T grows above the transition temperature, the chemical potential starts taking negative 
values and the fugacity becomes smaller than 1, thus moving the poles from the real axis 
into the complex plane and hence removing the previous singular behavior from all three 
responses. In this sense, the presence of these singularities in the non-condensate part of the 
response may be considered as an indirect signature of the existence of a Bose condensate. 

At fixed T, the low q response presents nonvanishing contributions coming from both 
the coherent and the incoherent parts of the response. Notice that, in contrast to what 
happens in the free Fermi gas caseB, the low g's coherent response may actually overcome 
the incoherent one. As q is risen, the coherent response decays rapidly to zero and thus the 
total response becomes mainly incoherent. As it is also apparent from the figure, this is 
almost the case when Q ^ 2, i.e., when q ?^ 2T^^'^. On the other hand, when q remains fixed 
and T is low, Q becomes large and hence S{Q, u) ^ SindQ, i^), that is, the total response 
becomes mainly incoherent. As the temperature is risen, coherent contributions become 



more relevant and eventually, depending on the value of g, at T :^ q both Sine and S, 



coh 



carry similar contributions to the total response. Therefore the evolution with T at fixed q 
can be understood as follows: at very low T, Q is large and so the response behaves as if it 
were in the high momentum transfer limit, where S{Q, v) ~ SindQ, v) and ScohiQ, v) -^ 0. 
As the temperature rises, Q shrinks and takes its lowest value at T = T^ before the Bose 
condensate disapears. When f ^ T^., the response still remains in the high momentum 



transfer limit, otherwise it enters into the low momentum regime. 

The behaviour of the three responses changes when the temperature rises above the 
Bose-Einstein transition temperature. In fact, at T > Tc the chemical potential is strictly 
negative and thus the fugacity z is smaller than 1. Under such circumstances, the poles of 
the responses that at T < 7^ were at z/ = itQ^ now move to v = iQ^ ± i2Q\/— Inz, and 
thus the singular points are removed from the real v axis, splitting each real pole into two 
complex ones. When either g or T increases, the distance from the poles to the real z/ axis 
grows and the shape of the response is smoothened. 

This behaviour is shown in Figs. (0), @ and (^ where the total, coherent and incoherent 
responses are depicted for several values of Q in comparison with the classical response 
Sci{Q, v'l T) computed from a Maxwellian momentum distribution 

SJQ, V- f) = J- e-^(^~'^)' , (12) 

^^ ' ^ 2QT3/2 0F ^ ^ 

which corresponds to the limiting case to which S{Q, u; T) should approximate when the 
temperature increases. 

As it can be seen from the figures, the way in which S{Q, z/; T) approaches SindQ, i^', T) 
(incoherent limit) and Sci{Q,i^',T) (classical limit) depends on the momentum transfer and 
on the temperature. At low T greater than T^, the total response approaches the incoher- 
ent one when q grows. However, both S{Q,v]T) and Sinc{Q,v;T) differ from Sci{Q,i^]T) 
because the momentum distribution at those temperatures noticeably depart from the gaus- 
sian shape of the classical n{k). When the temperature is risen, the difference between the 
momentum distribution in the classical and quantum quantum cases is reduced and so the 
high q response of the free Bose gas approaches the classical limit. Notice however that the 
difference between S{Q, z/; T) and Sci{Q, v, T) is still apparent at T = 20, which means that 
the classical limit is only well recovered at rather high temperatures. 

The evolution with T of the different responses can also be analyzed from the sum rules 
they satisfy. Sum rules are defined as the different energy-weighted moments of the responses 

^ticohi^) = / (^""Sinccohiq, i^; f) du 

J —oa 
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The first lower orders sum rules can be easily deduced and yield the following results 

m(°)(g;T) = no(T) [2n(g) + 1] +mW(g;r) = 5(g; f ) (14) 

^inll"?; ^) = no{T) + m^fl,nM T) = 1 (15) 

mS(g; t) = 2no(t) n(g) + mS^„,(g; t) = 5(g; t) - 1 (16) 

m(i) (g; t) = n,{f)f + m^ (g; t) = g^ (17) 

"^Snc(^; T) = no{f)f + m\l^c,nM T) = f (18) 

miii(g; T) = mlii,„,(g; T) = , (19) 

where ml^J{q) and ^(°„cco/i) nc('?) refer to the contribution to the sum rules coming from the 
integral of the non-condensate parts of the different responses. 

Eq. (|14|) yields the values of the static structure factor S{q] T) as obtained from the 
direct integration of Eq. (^). Two different contributions appear at temperatures below 
Tc. The first one results from the condensate term in S{q^ uj\ T) and so is proportional to 
noiT). The contribution of this term to the m^'^\q;T) sum rule grows as no{T) i.T/cf'j 
at low momentum transfer because at T < T^ the low q states close to the condensate 
are populated according to a T /(p law. Notice that this is an entirely coherent effect as 
can be seen from Eqs. (|1^) and (|16]). This behavour can be smeared out when particle 
interactions are allowed to take place as happens in liquid ^He at low temperatures. Also, 
as T~^/^5„c(g, ^] T) = SnciQ, ^', T) is a T-independent function of Q and z/ at T < Tc, the 
zero order i/-weighted moment of SndQ, ^) yields an universal curve from where both the 
evolution with q and the evolution with T of Snciq] T < Tc) can be extracted. According to 
the previous definitions M^^Q) = J z/"5(Q, u; f)dv = T"+3/2^W(g; y). 

M^J{Q) is depicted in the upper plot of Figure (||). At low Q, that is, at low momenta 
compared to T^/^, the divergences appearing in S{Q, z/) at z/ = iQ^ take most of the strength 
and as a result M^J{Q) increases. However, the singularities in S{Q, u) aX u = ±Q^ are of 



the logarithmic type and thus their //-weighted integrals of the response give always finite 
results consistent with Eqs. (|1^-(|I9|). On the other hand, in the large Q limit M^J{Q) 
approaches the constant value T'^/"^ ^ 0.058, a result that is consistent with m^^\q, T) going 
to 1 — uqIT) and the temperature dependence of the condensate fraction value reported in 
Eq. (^. When this value is reached the total S{q,T) equals 1. The high Q limit can be 
reached either rising q at fixed T or lowering T at fixed q. Therefore, no matter what the 
value of g is, at T — i> the static structure factor goes to 1, a fact that is consistent with 
S{q, u;T = 0) being completely incoherent and equal to 6 {uj — (p). When the temperature 
is risen, the total response receives contributions from the non-condensate terms and thus 
S{q;T) grows. Finally, when T reaches the transition temperature M^J{Q) equals M^^\Q) 
but its value depends on the the momentum transfer. If q is still high compared to T^^'^, 
the non-condensate part of the total response does not qualitatively differ from the T = 
case and thus the the static structure factor remains close to 1. Otherwise, the response 
enters into the low Q regime as seen in Fig. (|I]) and S{q; T) departs from unity. Taking into 
account that m^^\q;T) = 1 and that the coherent response of the free Bose gas is positive 
defined, the asymptotic value S{q; T) = 1 is always reached from above. 

When T exceeds the Bose-Einstein transition temperature the condensate fraction value 
decays to and the previous sum rules coincide with their non-condensate parts. The 
evolution with T of S{q] T) is sketched in the lower plot of Fig. @. 

As the logarithmic singularities in the dynamic structure function are smeared out at T > 
Tc, the limiting value S(g— ^0;T) now becomes finite. When q increases, S{q;T) decreases 
and asymptotically approaches 1 from above, thus indicating that the total response is 
reaching the incoherent limit in which S{q,oj;T) ^ Sinc{q,<^',T) and Scoh{(j,<^',T) ^ 0. In 
the high T limit S{q] T) ^ 1 at all g's and this reveals that the classical limit has almost been 
reached, as Sci{q',T) is known to be equal to 1 for all q. Moreover this limit is consistently 
reached when T -^ oo, as in fact it is straightforwardly seen from the very definition of T 
that this can be achieved either increasing the temperature at fixed density or reducing the 
density at fixed temperature. 



Another interesting feature of the high momentum transfer behaviour of the response is 
the scahng property, as it is well known that when q—>-oo the dynamic structure function 
does not depend anymore on q and uj separately but only through the West scaling variable 
y = {u/q-q)/M. 

In the free Bose gas case scaling may be reached in two different ways depending on the 
temperature. At T <Tc this can be done introducing a Compton profile J^dQ: ^j T) as 
follows 

J-(Q,r;T) = 2Q5(Q,z/;f) (20) 

where Y = {y/Q — Q) /2. With this definition, the non-condensate contribution to the 
incoherent response becomes temperature and momentum independent 

4c,nc(n = -^ln[l-e-^'] , (21) 

and so becomes an universal function valid for all values of Y as long as T < Tc. The non- 
condensate contribution to the total response reaches the incoherent limit at high momentum 
transfers, so JnciQi ^i T) approaches Jnc,incO^) when Q ^ c>o. This is indeed the case as 
shown in Fig. (|^). As can be seen, already aX Q = 2 the incoherent limit seems to be almost 
reached, which means that at any temperature below T^ scaling is well recovered when the 
transferred momentum is larger than q = 2yTc^5. 

Scaling can also be recovered at temperatures above the Bose-Einstein transition tem- 
perature. However in that case the following alternative definition of the Compton profile 
may be used 

J+iQ, Y; f) = 2Q f'/^S{Q, u- f) . (22) 

With this definition, the classical response becomes Q and T-independent according to the 
expression 

^ci(n = ^e-^\ (23) 



and therefore this is the hmiting case to which the response of the free Bose gas tends 
when the temperature is risen. The response J~^{Q,Y;T) and the incoherent response 
JticiQ^ F; T) at Q = 1 are compared at different T's with J'^i{Y) in Fig. (|^). Both functions 
approach the classical behaviour when the temperature is risen, even though it is remarkable 
how far should the temperature go in order to reach that limit. In any case and as expected, 
the incoherent response approaches faster J^(y), a fact that can be understood recalling 
that the total response contains also coherent contributions that do not scale in the same 
variables. In this way, the different velocity at which both functions approach J^iiX) gives 
a measure of how fast the coherent response goes to zero. Finally the main reason why the 
scaling in Eq. ( p2D is proportional to T^/'^ may be understood recalling that in the classical 
regime z is proportional to T^/^, so that in the end both dependencies cancel and J^iO^) 
becomes T-independent. This explains why the scaling law at T > T^ and T<Tc should be 
different, as in the latter case the fugacity vanishes and therefore it can not introduce any 
explicit dependence on T as happens at high temperatures. 

In summary, the dynamic structure function of the free Bose gas at finite tempera- 
tures is shown to contain both coherent and incoherent contributions. At low temperatures 
compared to the Bose-Einstein transition temperature T^ , the total response is mainly in- 
coherent and tends to the limiting incoherent value 6{uj — q'^/2m). When the temperature 
is risen, coherent contributions presenting logarithmic singularities at cj = ±g^/2r7i appear, 
and they actually overcome the contribution of the incoherent response when T ^ T^. When 
T is increased above T^, the logarithmic divergences disappear but a visible signature of their 
presence below Tc still remain in both Sine and Scoh- When T ^ T^ the coherent response 
cancels while the incoherent and total responses approach the classical limit. 

The incoherent limit is always recovered when the transferred momentum is high, while 
the classical limit is only reached at high T. Below T^, the presence of singularities in the 
response is also reflected in their lowest order energy-weighted sum rules. Finally, in the 
high q limit the total response can be written in the form of a Compton profile that scales 
in the West scaling variable Y . The temperature dependence can also be removed from 

10 



the response by using an approppriate redefinition of Y, even tliougli tlie transformation 
is different at T < Tc and at T > T^ due to tlie different temperature dependence of tlie 
cliemical potential. Despite tlie simplicity of the system analyzed and the fact that the 
excitation spectrum of an interacting bose system will be different the results described in 
this work are expected to enlight some aspects of the finite temperature response of both 
weakly and strong interacting Bose systems. 
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FIGURES 
FIG. 1. Non-condensate contributions to the total {S{Q, u), solid lines), coherent {ScohiQ, '^), 
dashed lines) and incoherent {Sinc{Q,i^), dotted lines), for Q = 0.7, 1.5 and 2, at any temperature 
below Tc. 

FIG. 2. Temperature dependence above Tc of the total {S{Q,i';T), solid line), coherent 
{ScohiQ I'^'jT), dashed lines) and incoherent {Sinc{Q,J^',T), dotted line), in all cases at Q = 0.7. 
The classical limit Sci{Q-,v',T) is also shown at each temperature (dot-dahsed lines). 

FIG. 3. Same plots as in Fig. (^) aX Q = 1.5. 

FIG. 4. Same plots in Fig. (1) at Q = 2. 

FIG. 5. Contribution of the non-condensate part of the response to the zero order sum rule. 
Upper plot: M^?(Q) at T < fc- Lower plot: m^^\q,f) at T > 7^, for T = 7 (solid line), T = 8 
(dotted line), T = 10 (short-dashed line), and T = 20 (long-dashed line). 

FIG. 6. Q dependence of the Compton profile JnciQ^ ^) at T < Tc compared with the universal 

inc,nc\ ) ' 

FIG. 7. Temperature dependence in the range T > Tc at fixed Q = 1 of total J'^{Q,Y;T) 
(left hand side) and the incoherent J^^{Q,Y;T) (right hand side), compared to the universal 
j;['(y) (solid line). 
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